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i Abstract 



q | We study the chiral symmetry in two-color QCD with Nf massless flavors at 

Q> 1 finite temperature, using an effective theory. For the gauge group SU(2), the 

chiral symmetry group is enlarged to SU(2iVy), which is then spontaneously 
broken to Sp(2A r j) at zero temperature. At finite temperature, and when the 
axial anomaly can be neglected, we find a first order phase transition occur- 
ring for two or more flavors. In the presence of instantons, the symmetry 
restoration unambiguously remains first order for three or more massless fla- 
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vors. These results could be relevant for lattice studies of chiral symmetry at 



^1" 

(N 



X 



finite temperature and density. 



ON : I. INTRODUCTION 

An important issue in QCD is the behavior of the theory under extreme conditions, 
Ph. such as a high temperature or large density. The high temperature region is relevant for 
an understanding of the early stages after the big bang, whereas high densities occur in 
e.g. neutron stars. In particular, with the upcoming events at RHIC such conditions may 
be obtained experimentally, thus providing a direct test of the theoretical predictions and 
increasing our knowledge about QCD at high temperatures and densities. 

Of special interest is the phase diagram of QCD, i.e. whether there are any phase 
transitions or cross-over regions when the temperature and chemical potential is increased 
from zero. For instance, QCD is well known to have a phase transition separating a zero 
temperature phase, where the chiral symmetry is broken, from a high temperature chirally 
symmetric phase. At zero chemical potential, the order of the finite temperature phase 
transition was studied in using the linear cx-model as an effective theory. Numerical 
studies on the lattice agree with the analytical results and indicate a second order phase 
transition for Nf = 2 massless flavors, which becomes first order for Nf = 3 massless or light 
flavors [[|. The best quantitative estimate for the critical temperature T c comes from the 
lattice simulations, giving a numerical value around T c ~ 150 MeV. 

Despite the time consuming algorithms for simulating full QCD on the lattice, the finite 
temperature observables are at least in principle calculable. However, when it comes to the 
introduction of a chemical potential, /x, the numerical simulations encounter severe problems. 
The principal difficulty is due to the functional measure (or more precisely, the determinant 
remaining after the fermions have been integrated out), which becomes complex-valued 
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for finite fi and thus makes the usual importance sampling impossible. Although several 
different proposals have been made in order to overcome this problem, such as the use 
Langevin algorithms || and the so called Glasgow method ||, the successes have so far 
been rather limited. A quantitative understanding of the chiral symmetry restoration in 
QCD at finite chemical potential is still lacking. 

Considering the above mentioned difficulties, it may very well prove worthwhile to study 
the finite chemical potential problem in a somewhat simpler context. In this respect, QCD 
with SU(2) replacing SU(3) as the gauge group is of considerable interest as a toy model. 
Specifically, since the four-dimensional, non-abelian structure is still maintained, some of the 
general properties in the simplified theory should be rather similar to the SU(3) theory. Thus, 
with such a minor modification one could hopefully obtain some insights of a qualitative 
nature into the structure of full QCD. 

However, although two-color QCD apparently resembles real QCD in many respects, 
there are nevertheless some striking differences. Obviously, hadron spectroscopy is going 
to be completely different. Another example is the consequence of a formation of diquark 
condensates, which has a different interpretation in two-color QCD and full QCD. In real 
QCD, a nonvanishing diquark condensate breaks the local gauge symmetry, and becomes 
an order parameter for color superconductivity ||, a new state in QCD which has received 
much attention lately (see e.g. 0] and references therein). In two-color QCD, the situation 
is slightly different since the diquarks may form a gauge singlet condensate. The condensate 
would still break baryon number conservation though, and hence become an order parameter 
for superfluidity instead. 

The main motivation for replacing SU(3) by SU(2) is the possibility to study the theory 
numerically. Since the fermions are in a pseudo-real representation for SU(2), the functional 
measure is real and positive, even with the introduction of a chemical potential. Therefore, 
two-color QCD at finite density can be studied on the lattice using standard methods. 
Preliminary lattice studies at finite chemical potential, both in the quenched approximation 
and with unquenched fermions ||, indicate a nonvanishing diquark condensate at any 
chemical potential \x > in the chiral limit. Due to the pseudoreal representation, when 
/i = the (qq) and the usual (qq) condensates can in fact be rotated into each other without 
any cost in energy. Moreover, analytical results using QCD inequalities and models 
based on instanton induced interactions |10[ agree with the overall picture emerging from 
the numerical results. 

There is thus an increasing amount of evidence in support of the above mentioned behav- 
ior of SU(2) QCD at finite //. However, very little is still known about the thermodynamics 
of two-color QCD, although some results were obtained rather recently O]. This paper is 
devoted to such a study of the chiral condensate at finite temperature, using an effective 
theory for the condensate, and the fluctuations around it, to describe the critical behavior. 
Our motivation is twofold; first of all, it is theoretically interesting to investigate the whole 
phase diagram of two-color QCD, for instance when it comes to the question of possible 
phase transitions and the connection to a superfluid phase. Secondly, it is of course also 
valuable to have a comparison between future numerical efforts in two-color QCD at finite 
temperature (and density), and other approaches based on different methods. Hopefully the 
combined methods will give a better understanding of two-color QCD, as well as a glimpse 
of the general features in full QCD. 
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The paper is organized as follows. In the next section we recall the chiral symmetry 
and the symmetry breaking pattern for two-color QCD, and write down an effective theory. 
Then in section III we use the effective theory to analyze the chiral symmetry breaking at 
finite temperature. We end in section IV by giving our conclusions and an outlook. 



II. THE CHIRAL SYMMETRY IN TWO COLOR QCD 

From now on, we assume that the gauge group is SU(2) and furthermore that there are 
Nf flavors of massless quarks. The part of the QCD action containing the fermions is given 
by, 

Nf N f 

s = E / M - wA;(m) 7 ^7 = -i E / **%mt . (!) 

/=1 /=1 

where the metric is g^ iv = diag(H ) and the color matrices T a = a a /2, with a a the usual 

Pauli matrices. Defining ay, = (1, a) and er M = (1,—a), we then split the four-component 
notation into the Weyl formalism as follows, 

Writing the action (|I|) in the Weyl two-component formalism we have, after a partial inte- 
gration, 

S=-iY,Ux [\ jf *» (d,5l + igA^T*)?) \) + fyf (d,5j - igA^T")?) £ j7 ] . (3) 
/=i 

Using now the pseudo-reality relation in SU(2), 

02<y a 02 = , (4) 

where a T denotes the transpose of a, the interaction term for A can be written as, 

\ jf (T a )l\) = A/TjA/ = -(X f a 2 )T a (a 2 X f ) . (5) 

Since A/ transforms as 2 under a gauge transformation, 02 A/ transforms as £, i.e. like 2. 
Thus, by defining 

n J £/ f 1 -V r> f c < / = l, . . . ,N/ 

^/ \a 2 X f / = N / + 1,...,2N / ' V/ \ A/a 2 / = N / + 1,...,2N / ' W 
the action in equation (|3|) finally becomes, 

S = E d 4 xQ f (d,- igAp*) a»Q f . (7) 
/=i 

Hence, the naive symmetry SU(iV/)iXSU(iV/)flXU(l)B is enlarged to an SU(2iVy) symme- 
try, which becomes manifest when the action is rewritten as in equation (|7|). It should be 
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noted that the U(l)s symmetry becomes a subgroup of SU(2iV/), and also that the deriva- 
tion only holds true at a zero chemical potential /i, since at finite /i the SU(2iVj) symmetry 
is explicitly broken to the original 811(^)^x811(^)^x11(1)^ symmetry |5). 

Now consider a spontaneous breaking of the SU(2iVy) chiral symmetry [12]. Restricting 
ourselves to a Lorentz and gauge invariant condensate bilinear in Q, the structure has to be, 
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where q are the color indices. In order to obtain overall antisymmetry the flavor structure 
must also be antisymmetric, and hence the chiral condensate S must be of the form, 



S OC Qi Cia Qj C2 /3€ 

where E is a (2Nf) x (2Nf) matrix, 



1 
-1 



(9) 



(10) 



with and 1 Nf x Nf matrices. If X ^ 0, the chiral symmetry is spontaneously broken and 
E is only invariant under a subgroup of SU(2iVy), specified by the requirement that for any 
element H in this subgroup, 



Q T H T EHQ = Q T EQ . 



(11) 



The set of matrices H satisfying equation (|TT|) form the symplectic group Sp(2iV/), so for 
a non-zero value of S, the symmetry breaking pattern is 3U(2Nf) — > Sp(2A r /). Since 
Sp(2Nf) contains (2Nj + Nf) independent parameters, the number of Goldstone bosons are 



(AN] 



1) - (2N 2 f + N 



f) 



2Nf 



N f -1. 



To describe the symmetry breaking and the corresponding low energy degrees of freedom 
we only need to take into account the condensate and the massless fluctuations around it, i.e. 
the Goldstone bosons. The coset space is given by SU(2iV/) /Sp(2iV/) and the effective action 
can be parametrized by a (2Nf)x(2Nf) complex antisymmetric matrix $, with (ANj — 2Nf) 
independent real-valued parameters, 



(12) 



Under an SU(2iVj) transformation M, $ transforms as $ — > M$M T , and the antisymmetry 
of $ is obviously maintained. 

The most general, renormalizable effective Lagrangian for invariant under the SXJ(2Nf) 
symmetry is, 
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A 2 Tr 
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(13) 



If one adds also a vector field, it is possible to enhance the global symmetry SU(2iV/) to 
SU(2iV/) xSU(2iVy), a fact that was explored in [13|. In this article we will, however, only 



consider the chiral symmetry breaking and the associated Lagrangian given above. 

Clearly, the Lagrangian in equation ( |i~3D is invariant under U(2Nf), and not just 
SU(2iVj). Additional terms that only respect the SU(2iVj) symmetry appear when the 
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axial anomaly is taken into account. The corresponding lowest-dimensional term for an 
even, antisymmetric 2n x 2n matrix Q is the Pfaffian, Pf(f2) ||14|| , given by 



pf (n) = J2(-i) p n ili2 n i3ii ■ ■ ■ n i2n _ li2n (i4) 

z n. p 

where z& = 1, . . . , 2n, and P denotes a summation over all permutations of {«!,..., i2 n }, 
with (-l) p the sign of the permutation. The square of Pf (O) is simply the determinant of 

a 

The action invariant only under SXJ(2Nf) finally becomes, 

S = J d*x {L -c [Pf ($) + Pf($ t )] } = J d 4 x {L + L P f.} , (15) 

where, in four dimensions, the different coupling constants have mass dimensions [Aj] = 0, 
[to] = 1 and [c] = 4 - N f . 

To study the symmetry breaking pattern SU(2iVj) — > Sp(2iV/), <£> should acquire a non- 
zero expectation value $o of the form 

^o = 0o(_° 1 I) , (16) 

where 1 is the Nf x Nf unit matrix and 0o a constant. The matrix structure in equation 
(|i~6|) clearly respects the Sp(2Afy) symmetry, and when m 2 < the potential in the action 
(p"5|) is minimized for <po ^ 0. Furthermore, stability at large values of $ requires A2 > 
and Ai + \2/2Nf > 0. With this symmetry breaking term it is straightforward, although 
rather tedious, to check explicitly that there are exactly 2Nj — Nf — 1 (2Nj — Nf) number 
of Goldstone bosons for c 7^ (c = 0). 



III. CRITICAL BEHAVIOR AT FINITE TEMPERATURE 

In this section we will study the influence of a finite temperature T on the chiral conden- 
sate, and the order of the phase transition to the SU(2Nf) symmetric phase as T reaches the 
critical temperature, T — > T c . At zero temperature, m 2 < and we have 0o 7^ 0, breaking 
the SU(2iVf) symmetry spontaneously to Sp(2iVj), as stated in the previous section. We 
will first assume that c = 0, so that the chiral symmetry is enlarged to U(2iVj), and later 
discuss possible consequences of a non-zero c. 

The critical behavior at finite temperature is determined by the infrared fixed points 
in the coupling constant space of Ax, A 2 and the mass parameter. Hence, to study the 
renormalization-group flow we need to calculate the scale dependence for the these parame- 
ters. Although it is of course possible to obtain the /3-functions from a direct calculation of 
the relevant Feynman diagrams, we will here use the simpler method of the effective action. 
Needless to say, the two different methods are completely equivalent. 

When the /3-functions are calculated from the effective action they will correspond to 
a system in (4 — e) dimensions, where e C 1. But since the extension of the timelike 
direction is only (3 = l/T (in the imaginary time formalism), the large scale fluctuations in 
the vicinity of the critical temperature is governed by an effective action in three dimensions. 
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Nevertheless, the e-expansion has proved to be a useful tool to describe critical behavior in 
three dimensions, where e = 1, both for the chiral symmetry restoration in real QCD Q but 
also, and to a great numerical accuracy, in several different condensed matter systems (see 
e.g. 
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In the effective action, we need to distinguish between the interaction terms with cou- 
pling strength Ai and A2 respectively. To accomplish that, the constant <f>o will be labeled 
differently for different matrix elements of ($0)^ 

O i = 1, . . . , (N f - 1) , ./ (A/- • 1) (2N f - 1) . 

foi = N f ,j = 2N f . 

-<i> a i = (N f + l),...,(2N f -l), j 1 (A/ 1). { <) 

-0 6 i = 2N f ,j = N f . 

Using equation (|TT[), assuming that (fi Q is real-valued and expanding the action fllS] ) to 
one-loop order around $0 we find, 



Si-ioop = / d 4 x |-ip T (d 2 + m 2 ) p - l -p T M int .p - m 2 [(N f - 1)0 2 + 



,2 , ^ 



- 4Ax (N f - 1)# + 2 - 2A 2 (JV> - l)0 4 a + 



2 



,4 , ^ 



= J d*x [-\p T (9 2 + m 2 ) p - l -p T M- int .p - V ] . (18) 

In the expansion above, we have put all the Nf(2Nf — l) fluctuations of the complex matrix $ 
into a real- valued, 2Nf(2Nf — l)-component vector p, where Mj nt . contains all the nontrivial 
terms depending on Ai and A2. In the last line, Vo is nothing but the classical potential, and 
by disregarding the different labelings and replacing (f) a = <Pb = <Po, h has a minimum for 



* = wU) ( - m2>0) - (19) 

Since M int can be made symmetric, it can be diagonalized by an orthonormal matrix P, 
M D = PM int P T = diag(m 2 , m 2 , . . . , ml Nf{mf _ 1} ) , (20) 

and by changing variables to p = Pp and m 2 = m\ + m 2 , k = 1, . . . , 2Nf(2Nf — 1), we find 
the following Euclidian effective potential after a gaussian integration, 

V eS = Vo + 7^7 E LogDet(-9| + m 2 ) , (21) 
21/4 fc 

where V4 is the four-dimensional volume and the subscript E denotes a Euclidean metric 
V = diag(+ + ++). 

Now, as the temperature T is turned on there is a thermal correction to the zero 
temperature mass parameter m, and when the temperature reaches the critical temperature 
T c we have 0o — > 0, or in other words m 2 +m 2 - c — > 0. Hence, the critical behavior occurs when 
the total mass m 2 + m\ is fine tuned to zero, as usual in critical phenomena. Therefore, in 
the vicinity of the fixed point we can ignore the m dependence in the effective potential ( pT|) . 
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and keep only the remaining divergent pieces in the second term to calculate the /^-functions 
for Ai and A2. We find, using dimensional regularization in (4 — e) dimensions, 



V * = V >-5ftZ«A. (22) 

k 

Finding all the eigenvalues of M int ., for arbitrary Nf, seems a formidable task. Fortu- 
nately, this is not necessary, since 

J2 m k = TrMf, = TrM? nt . (23) 

k 

Thus, we can use the matrix Mj nt . directly in the effective potential, 

V eS = V - ^TrM^, . (24) 
Finally, by writing down the explicit form of M int . the /3-functions become, 
A = ^ = -eA, + (^1) A? + A,A 2 + (£) A 

ft=f^ =-^ + (%5)A|+(4 r )A 1 A 2 , 



2 

2 J 



(25) 



where \i is an arbitrary mass scale. 

As is well known, an infrared fixed point A* corresponds to $ = and positive eigenvalues 
of the matrix Wij = (d/3i/d\j)\\=\*, for i, j = 1, 2. From equation (p5|) we find that the only 
infrared stable point to order e is given by 



A ^ 2j V|-7V / + 4 ' A ^°' < 26) 



and only for N f < (1 + V17) /4 = iV^ ax ~ 1.28. When JV) > A^f ax the fixed point is infrared 
unstable along the A2-axis. 

Assuming that the predictions from the e-expansion remain valid in three dimensions 
and in addition that the theory evolves from a point within the stability region (A2 > 0, 
Ai + \2/(2Nf) > 0), we conclude that there is a second order phase transition in two-color 
QCD with Nf < (1 + vTT)/4 massless flavors, to a phase where the spontaneously broken 
XJ(2Nf) symmetry is restored and the chiral condensate disappears. Physically, the only 
relevant case is Nf = 1 where there is just one independent, complex element in $. The 
phase transition will thus be characterized by the critical exponents of the 0(2) linear o- 
model. For two or more massless flavors, the instability in the direction of A2 should result in 



a first order phase transition induced by fluctuations [16] , the analogy in (4 — e) dimensions 
of the Coleman- Weinberg mechanism [[H] . 

Let us now consider the effects of having a non-zero value for c around the critical 
temperature. In that case, we have to add to the Lagrangian the term 



L Pf . = -c Pf ($) + Pf ($ T ) . (27) 
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For a single flavor this is a linear term, so c acts like an external magnetic field. As a result, 
the symmetry is never really restored, since for c ^ the potential always favors a non- 
zero expectation value. For two flavors, the symmetry group is SU(4)~0(6) and only half 
the number of fields become massless at the transition point. The phase transition could 
therefore be of second order and fall into the 0(6) universality class. 

Adding a third flavor, the term proportional to c becomes cubic, which in itself drives 
the phase transition first order. For Nf > 4, the operator is either marginal (Nf = 4) 
or irrelevant (Nf > 4) in four dimensions. In both cases, the critical behavior is already 
determined by Ai and A 2 , so the transition remains first order, induced by fluctuations. 



IV. CONCLUSIONS 

We have presented a study of the chiral condensate in two-color QCD at finite temper- 
ature and zero chemical potential. Based on an effective theory, we find conclusive evidence 
for a high temperature region where the condensate disappears. Without the anomaly, there 
is a second order phase transition for one flavor, and a first order phase transition induced 
by fluctuations when Nf > (1 + VT7)/4. 

When the anomaly is taken into account, the situation becomes rather different for 
Nf < 2, since there is strictly speaking no phase transition for Nf = 1 and possibly a second 
order phase transition for Nf = 2. When Nf > 3, the phase transition is still first order, 
but for Nf = 3 it is not induced by fluctuations. 

The order of the phase transition for a small number of flavors is therefore quite sensitive 
to the value of the instanton-induced term c, a situation rather similar to the one found for 
the chiral symmetry in full QCD 

As with three colors, instantons in two-color QCD become less significant as the tem- 



perature is increased, since the instanton size is cut off by the temperature [f[q| . Hence, the 
temperature dependence of c is what really matters. At T = the anomaly term should 
of course be present, but it could happen that c decreases with an increasing temperature, 
and in such a way that it does not play any significant role at the temperature scale where 
the condensate disappears. In that case, the critical behavior is governed by a restoration 
to the XJ(2Nf) symmetric phase. The only exception occurs for Nf = 1, where any relic 
of the linear term destroys the second order phase transition at c = 0. If, on the other 
hand, c is not neglectable, the symmetric phase corresponds to a restoration of the SU(2iVj) 
symmetry. 

Keeping in mind that we have considered a simplified version of QCD, it is somewhat 
reassuring to find a behavior reminiscent of the chiral symmetry in full QCD at finite T, in 
particular since both the original symmetry and the remaining subgroup are rather different 
in the two cases. This may indicate that some of the qualitative features of two-color QCD 
can indeed be taken over to real QCD, at least with some care. 

Another interesting question is how the results derived here are related to the finite 
temperature behavior of two-color QCD at \i > 0. At zero temperature, the formation 
of a diquark condensate seems to be well established by now, and in the chiral limit the 
condensate is nonvanishing at all chemical potentials, which implies that the U(l)s part of 
the overall SXJ(Nf) L xSXJ(Nf) R xXJ(l) B symmetry at finite fi is spontaneously broken (the 
symmetry breaking pattern is SXJ(N f ) L xSXJ(N f ) R xXJ(l) B ->Sp(iV / )xSp(iV / ) 0). Recent 
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lattice simulations at nonzero quark masses m |19[ seem to indicate that the chiral condensate 
is nonzero at fi = and then decreases smoothly as the chemical potential becomes larger, 
whereas the diquark condensate is zero up to some /x and then increases. Most likely, there is 
a genuine phase transition separating the two regions at some critical chemical potential /x c , 
and there is presumably also a mixed phase for /io ~ ^ ~ /^c, where none of the condensates 
vanish. 

Thus, even though the SU(2iV/) is explicitly broken by both a finite mass m and a 
finite chemical potential /a, the chiral condensate is still nonvanishing for small values of m 
and ft, as long as ll < \i c ~ \fm. Consequently, if both m and ll can be treated as small 
perturbations, it is not inconceivable to assume that the chiral phase transition at finite 
T can be modified to also include the case m, ll > 0. For example, a second order phase 
transition could be changed into a cross-over region, whereas a first order one presumably 
would become weaker, but remain. By continuity arguments, the chiral phase transition (or 
crossover) could also meet the phase transition separating the chiral and diquark condensate, 
possibly at a critical or tricritical point. Furthermore, since there is no condensate at all 
above the finite T phase transition presented in this letter (at ll = 0), this could indicate 
that there is yet another phase transition at finite T and ll, between the superfluid and 
non-superfluid phases. 

Although the above comments are consistent with the lattice results obtained so far, it 
should be emphasized that the remarks are highly speculative at the present stage. Hopefully, 
the combined results of future analytical and lattice investigations will be able to clarify the 
situation. 
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